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We analyze the Brownian thermal noise of a multi-layer dielectric coating, used in high-precision optical 
measurements including interferometric gravitational-wave detectors. We assume the coating material to be 
isotropic, and therefore study thermal noises arising from shear and bulk losses of the coating materials. We 
show that coating noise arises not only from layer thickness fluctuations, but also from fluctuations of the 
interface between the coating and substrate, driven by internal fluctuating stresses of the coating. In addition, 
the non-zero photoeleastic coefficients of the thin films modifies the influence of the thermal noise on the laser 
field. The thickness fluctuations of different layers are statistically independent, however, there exists a finite 
coherence between layers and the substrate-coating interface. Taking into account uncertainties in material 
parameters, we show that significant uncertainties still exist in estimating coating Brownian noise. 

PACS numbers: 



I. INTRODUCTION 

Brownian thermal noise in the dielectric coatings of mir- 
rors limits some high precision experiments which use optical 
metrology. This thermal noise is currently a limit for fixed 
spacer Fabry -Perots used in optical clock experiments [1] and 
is estimated to be the dominant noise source in the most sen- 
sitive band of modern gravitational wave detectors (e.g., ad- 
vanced LIGO, GEO, advanced VIRGO andLCGT) [2-6]. Re- 
cent work has indicated the possibility of suppressing the var- 
ious kinds of thermal noise by redesigning the shape of the 
substrate and the structure of the multi-layer coating [7, 8]. In 
this paper, we seek a comprehensive understanding of coating 
Brownian noise, by identifying individual sources of fluctua- 
tions, calculating their cross spectra using the fluctuation dis- 
sipation theorem [9, 10, 13], and finally evaluating how each 
of the sources and their correlations add up to the total noise. 

As a starting point, we will assume each coating layer to be 
isotropic, and hence completely characterized by its complex 
bulk modulus K and shear modulus — each with small imag- 
inary parts related to energy loss in bulk and shear motions. 
The complex arguments of these moduli are often referred to 
as loss angles. While values of K and jj. are generally known, 
loss angles vary significantly according to the details of how 
coating materials are applied onto the substrate and their com- 
position. Since the loss angles are small, we will still use K 
and ji to denote the real parts of the bulk and shear moduli, 
and write the complex bulk and shear moduli, K and p. as 



K = K{l+i<j)B), A=M(1+%)- 



(1) 



Here we have used subscripts B and S to denote bulk and 
shear, because these will be symbols for bulk strain and shear 
strain. Note that this definition differs from previous literature 
and measurements, which used 0|| and to denote losses in- 
duced by elastic deformations parallel and orthogonal to the 
coating-substrate interface [11]. We argue in Appendix C that 
0|| and 0x cannot be consistently used as the loss angles of a 
material. 

Brownian thermal fluctuations of a multilayer coating can 
be divided as follows: (i) thickness fluctuation of the coating 
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Figure 1 : Schematic plot of a mirror coated with multiple dielec- 
tric layers. Shown here are the various fluctuations that contribute to 
coating noise, i.e., fluctuations in the amplitude and phase of the re- 
turning light caused by fluctuations in the geometry (e.g. layer thick- 
ness Slj, layer area stretch {8A/A)j and interface height Zs) of the 
coating-substrate configuration and in refractive indices [dnj{x,y,z)] 
of the layers. 



layers, (ii) fluctuation of the coating-substrate interface, and 
(iii) refractive index fluctuations of the coating layers associ- 
ated with longitudinal (thickness) and transverse (area) elastic 
deformations — as illustrated in Figure 1 . Using what is some- 
times refeiTed to as the Levin's direct approach [10] (based 
on the fluctuation dissipation theorem), one can lump all three 
contributions into one quantity, and calculate its noise spectral 



2 



density by calculating the mechanical dissipation rate when a 
distribution of mechanical forces are applied at various loca- 
tions on the coated mirror [as has been done by Vyatchanin 
et al. [12]]. However, in order to obtain insights into coat- 
ing noise that have proven useful we have chosen to calculate 
the cross spectral densities for each of (i), (ii), and (iii), and 
provide intuitive interpretations of each. We will show, that 
(i) and (ii) above are driven by both bulk and shear fluctua- 
tions in the coating, in such a way that thickness fluctuations 
of the j-th layer 5lj, or in transverse locations separated by 
more than coating thickness, are mutually statistically inde- 
pendent, yet each 5lj is correlated with the fluctuation of the 
coating-substrate interface Zs — because Zs is driven by the sum 
of thermal stresses in the coating layers. We will also show 
that when coating thickness is much less than the beam spot 
size, the only significant contribution to (iii) arises from lon- 
gitudinal (thickness) fluctuations. 

This paper is organized as follows. In Sec. II, we express 
the amplitude and phase of the out-going field in terms of 
fluctuations in the coating structure, thereby identifying the 
various components of coating thermal noise. In Sec. Ill, we 
introduce the loss angles of isotropic coating materials, and 
use the Fluctuation-Dissipation Theorem to calculate the cross 
spectral densities of the coating thermal noise ignoring light 
penetration into the multi-layer coating. In Sec. IV, we dis- 
cuss in detail the cross spectra of all the components of the 
coating structure fluctuation, thereby obtaining the full for- 
mula for coating thermal noise, taking light penetration into 
account. In Sec. V, we discuss the effect of light penetration 
on coating thermal noise, using typical optical coating struc- 
tures. In Sec. VI, we discuss the dependence of thermal noise 
on the material parameters, and optimize the coating structure 
in order to lower the thermal noise. In Sec. VII, we discuss 
how only one combination of the two loss angles have been 
measured in past experiments, and how other different combi- 
nations can be measured using a new experimental geometry. 
Finally, we summarize our main conclusions in Sec. VIII. 



II. COMPONENTS OF THE COATING THERMAL NOISE 

In this section, we express the coating thermal noise in 
terms of the elastic deformations of the coated substrate. 



A. Complex Reflectivity 

As illustrated in Figure 1, we consider a laser field nor- 
mally incident (along the —z direction) onto the mirror, with 
complex amplitude profile Uin{x,y) at a fixed reference plane 
(dashed line in the figure) and intensity profile I{x,y) = 
\uin{x,y)\'^. Henceforth in the paper, we shall use arrows (e.g., 
x) to denote the 2-dimensional vector {x,y) in the transverse 
plane, and boldface letters (e.g., x) to denote 3-dimensional 
vectors. 

Because coating thickness is much less than the beam spot 
size, the reflected field (traveling along the +z direction) at 



transverse location x has an amplitude given by 



Mout(^) =Ptot(^)Min(^), 



(2) 



which only depends on the complex reflectivity Ptot(^) and 
the complex amplitude of the incident field Uia{x), at the same 
location x — assuming no incident light from the substrate 
(i.e., 52 = 0). Here Ptot(^) can be separated into three factors, 



Ptot{x) = 



Vlix) 



Uin(x) 



vi{x) 



(3) 



in which vi{x) is the incident complex amplitude at the 
coating-air interface, while V2 (x) is the reflected complex am- 
plitude at that interface. 

The first two phase factors on the right-hand side of Eq. (3) 
are gained by the light when traveling across the gap between 
the fixed reference plane (see Fig. 1) and the coating-air in- 
terface; we therefore obtain, up to a constant phase factor. 



"out(-^) 

V2ix) 



Min(x) 



(4) 



where ko = (Oo/c is the wave number of the laser (cOq its angu- 
lar frequency) in vacuum, Zs{x) is the vertical displacement of 
the coating-substrate interface (from its zero point), and lj{x) 
is the thickness of the j-th coating layer — both evaluated at 
a transverse location x. 

The remaining complex reflectivity V2{x)/vi (x) can be de- 
termined as a function of the phase shift experienced by the 
field in each layer, as well as the reflectivity of each interface, 
as described in detail in Sec. V. We can write: 

vi/vi ^ p[(j)i{x),...,(j)N{x);roi{x},...,rN>;{x)] (5) 

Here p is the complex reflectivity of a multi-layer coating, 
measured at the coating-air interface, which in turn depends 
on the optical thickness (l)j{x) of each layer (j = 1,...,N) 
and the reflectivity rp p+i(^) = rp{x) of each interface, (p — 
0,. ..,N, with p ~N+1 representing the substrate, and p = 
the vacuum outside the coating). Assembling the above equa- 
tions (3)-(5), we obtain: 



0:o:{x) ^ e-"'° [^-«+^-' p [{0,(x)}; )}] 



(6) 



Brownian thermal forces lead to fluctuations in both the real 
and imaginary parts of this complex reflectivity. Fluctuations 
in the argument of the complex reflectivity phase modulates 
the out-going light and directly produces sensing noise. Fluc- 
tuations in the magnitude, on the other hand, amplitude mod- 
ulate the out-going light, and produces a ponderomotive force 
noise. 



B. Thermal Phase and Amplitude Noise 

Brownian thermal fluctuations in coating geometry and re- 
fractive index modify the total reflectivity Ptot(^) defined in 
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Eq. (6). The real and imaginary parts of 

Sptotix) 



5logptot(x) 



Ptot(^) 



(7) 



encode the amplitude/intensity and phase fluctuations of the 
reflected Hght at position x on the mirror surface. In particular, 
intensity fluctuation of the reflected light is given by 

m^2'lB^=2Re[5lo,p..,m (8) 
I{x) \Ptot{x)\ 

while phase fluctuation is given by 

5(l){x) = 5 arg [pm{x)] = Im [5logptot(x)] (9) 

In this way, if we further write 

^(x)-/C(x) = --J^5[logptot], (10) 

with both £, and ^ real-valued functions of x, with the di- 
mensionality of displacement; they will represent phase and 
amplitude noise, respectively. In particular, from Eq. (9), we 
have 



2ko^{x) = 8^{x), 



(11) 



and this means ^{x) corresponds to the spurious displace- 
ment measured by the reflected light due to phase fluctuations 
caused by the coating. 

The quantity ^ is connected to amplitude/intensity noise via 



2koC{x) = Re [Slogptot] = 



5l{x) 
21{x) 



(12) 



As we shall discuss in Sec. HE, C, will cause a fluctuating 
force on the mirror, and contribute to measurement noise, al- 
though the effect will be small for gravitational-wave detec- 
tors. 

Inserting the dependence of Ptot on p, Ij and Zs, we obtain 



m-i^x) = 



-8z.{x)-Y^5l,{x) 

1=1 

N 



dlogp 
"<9logp 



5(1) j{x) 



5rp{x 



(13) 



The first two terms are due to the motion of the coating-air 
interface at location x and thickness fluctuations of the lay- 
ers, while the last two terms are due to light penetration into 
the coating layers. In particular, the third term is due to fluc- 
tuations in the total phase the light gains when propagating 
within the y-th layer, while the fourth term is due to the (ef- 
fective) reflectivity of the p-th interface (with p = Q indicating 
the coating-air interface), whose origin will be explained be- 
low. 



C. Fluctuations 8(^j and 8rp 

Light propagating within the coating layers are affected by 
the photoelastic effect, namely an isothermal fluctuation in 
5nj{x) (note here that x is a 3-D vector) due to fluctuating 
Brownian stresses exerted onto the coating materials. Assum- 
ing isotropy of the coating materials, we can write 



5«,(x) = /3j-5,,(x)+j3/[5„(x)+5,,(x)] (14) 



with 



(9 log/ 



i3 - 



drij 
5 log A 



(15) 



Here L stands for longitudinal, and T stands for transverse, 
and the subscript Aj and Ij indicate fixing transverse area and 
longitudinal length, respectively. We have also used the usual 
strain definition 



1 



dxj dxi 



(16) 



where m,(x), i — 1,2,3 are components of the displacement 
vector of the mass element at position x. Please refer to Ap- 
pendix B for more details in defining the elasticity quantities, 
and Appendix A 1 for more details on the photo elastic effect. 

We note that in Eq. (14) S„ is the fractional increase in 
length (i.e., linear expansion) in the longitudinal direction, 
while Sxx + 5'vy is the fractional increase in the transverse area. 
According to Appendix A 4, we can ignore the second term 
representing area fluctuations in Eq. (14) when the beam spot 
size is much larger than the coating thickness. In this case, 
we write jSy in place for ^j, whose value can be expressed 
in terms of a particular component of the photo elastic tensor, 
see Eq. (A5). 

As we discuss in Appendix A 2, the (surviving) first term of 
Eq. (14) causes two effects for light propagating along each 
direction (i.e., +z and — z): it adds an additional phase shift 
onto the light, and it back-scatters a fraction of the light into 
the opposite direction. As we show in Appendix A3 [c.f. 
Eqs. (A16)-(A18)], theses effects can be described by mod- 
ifying the phase shift of each coating layer and the reflec- 
tivity 5rj of interface: 



5rj = kot]lij5l}. 
Here we have defined 



1-r 



2r 

1 +'•./-! 



2r 



dfj = - / S„{^j+\+z)co?,{2kQnjz)dz 



5Pj = -y S.^{zj+i+z)sin{2konjz)dz 



(17) 
(18) 

(19) 
(20) 
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for j> 1, 5/g = S/g =0, and 



TV 



(21) 



marks the z-coordinate of the top surface of the y'-th layer We 
can also write 



51 i 



Note that 



total coating 
thickness 



Szz{zj+\+z)dz. (22) 



= Zl >Z2 > ••• >ZAf+l =0 (23) 



Note that 5rj, as well as the last two terms in are due 
to back-scattering, and have not been considered by previous 
authors. 

Inserting Eqs. (17), (18) into Eq. (13), we obtain: 



m-im=-zs{x)-i^ 



f" 









u„{x,z)dz 
(24) 



where 



l-rjd\ogp 
l+rjdlogp 



2rj d(j) 



,7+1 



cos[2kfinj{z — Zj)] 



t]l5j^^^sm[2konj{z-Zj 



(25) 



a term that accounts for all effects associated with light pene- 
tration. Here we need to formally define 



5logp 



= 



(26) 



w+i 



since (pN+i does not really exist. Alternatively, we can also 
write formulas separately for £, and ^, using only real-valued 
quantities. For we have, 

m = -Zsix) 



- £ 51 J ix) + 3rf'5l'j (x) + ^^'51] {x 



where 



3": 



4 V d^j 



4 V 

Pjtl^^fdlogp 




,(27) 



(28) 



(29) 
(30) 



are transfer functions from the various 5/'s to the 
displacement-equivalent thermal noise. For ^, we have 



where 



3f5ljix) + .!7^'5l)[x) + ^'5l^ix) 



3'- 



2 

4 V ^^1 

4 V 

M^jdiogp 




(31) 



(32) 



(33) 
(34) 



Although for an arbitrary stack of dielectrics, ^ is compa- 
rable to the part of £, [c.f. Eq. (25)] that involves light pene- 
tration into the layers. In practice, for highly reflective stacks, 
the real parts of dlogp/d(j)j and dlog p/drj all turn out to 
be small, and therefore fluctuations in ^ should be much less 
than fluctuations in ^ . 



D. Mode selection for phase noise 

So far we have calculated phase and amplitude noise as 
functions of location x on the mirror surface. However, there 
is only one displacement noise that the light will sense. In 
this and the next subsection, show how ^{x) and ^(x) should 
be converted into measurement noise. In doing so, we rec- 
ognize that only one spatial optical mode is injected on reso- 
nance in the optical cavity, and this mode has a complex am- 
plitude of uo{x) at the mirror surface. Now suppose we have 
"in = uq{x) incident on the mirror surface, we will then have 
Uout{x) = Ptot{x)uo{x), which contains not only the resonant 
mode, but also other modes, which do not resonate in the cav- 
ity. 

Let us select only the component of Mout(^) resonates, then 
we have a complex reflectivity of 



P = 



[ u^-,{x)uout{x)d^x _ J pioi{x)I{x)d^x 



J u^uodx 



JI{x)d^x 



(35) 



specifically for the resonant mode, and hence independent of 
X. Here we have defined I{x) = \u(){x)\-^. Note that the bar on 
top of p represents averaging over the phase front, instead of 
averaging over time. 

Now, inserting Eq. (10) as definitions for ^ (x) and ^{x) into 
Eq. (35), we obtain the fluctuating part of p 



dp 
P 



2ih){t,~iQ, 



where 



J^{x)I{x)d^x 
JI{x)d^x ' 



JC{x)I{x)d^x 
JI{x)d^x 



(36) 



(37) 
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Note that likoE, is the additional phase gained by the return- 
ing Hght, while 2ko^ is the relative change in amplitude [see 
discussions in Sec. II B]. Focusing first on |, we note that this 
creates the same phase as that gained by the reflected light 
if the mirror does not deform but instead moves by ^ . In this 
way, ^ is an error in our measurement of the mirror's displace- 
ment. 



Nevertheless, as we will show in Sec. VB, the minor ampli- 
fication factor here is not enough in making amplitude noise 
significant, because ^ tend to be much less than ^, for the 
coatings we consider 



III. THERMAL NOISE ASSUMING NO LIGHT 
PENETRATION INTO THE COATING 



E. Conversion of Amplitude Noise into Displacement 

The amplitude thermal noise produces spurious GW sig- 
nal by modulating the radiation pressure acting on the mirror, 
which in turn drives spurious mirror motion. Let us first con- 
sider a single-bounce scenario, in which an incoming beam 
with intensity profile I{x), unaffected by thermal noise, is re- 
flected with an intensity profile I{x) + 8l{x), with 5l{x) in- 
duced by amplitude thermal noise. In this case, the mirror 
feels a thermal-noise-induced recoil force of 



^single 



'J^,^x. 



Using Eqs. (12) and (37), we obtain 



psingle _ 4/o^o - 



(38) 



(39) 



with /o the power incident on the mirror. If the mirror is 
within a cavity, then we need to consider both the increase 
in the circulating power (which we denote by 7^) with respect 
to the input power, and the coherent build-up of amplitude 
modulation within the cavity. We also note that now both the 
incident and reflected beam contains amplitude modulation, 
and that we must also consider the effect of this amplitude 
modulation on the input mirror. 

If we restrict ourselves to a single optical cavity on reso- 
nance, then the force thermal noise below the cavity band- 
width is given by 



^th 



16^0/. 



(40) 



Here is the circulating power in the arm cavity. Suppose 
both input and end mirrors have the same mass M, then the 
spectrum of cavity length modulation driven by the amplitude 
thermal noise is given by 



^amp 



32a)o/, 



(41) 



Note that ^ has the units of displacement, and therefore the 
pre-factor in front of in Eq. (41) is a dimensionless con- 
version factor from ^ to displacement noise. For Advanced 
LIGO, this cannot be completely dismissed at this stage, be- 
cause 



320)0/, ^ Ic 40 kg 

iQ.^c^^/Ti ' 800 kW ' m 



lOHz 

D./{27t) 



0.03 



(42) 



In this section, we compute the coating Brownian noise as- 
suming that the incident light does not penetrate into the coat- 
ing. This means light is promptly reflected at the coating- 
air interface, and therefore we should only keep the first two 
terms on the right-hand side of Eq. (13), which leads to = 0. 
We therefore consider only coating phase noise £, , in partic- 
ular its weights average throughout the mirror surface, £, , see 
Eq. (37). 



A. The Fluctuation-Dissipation Theorem 

The Fluctuation-Dissipation Theorem relates the near- 
equilibrium thermal noise spectrum of a generalized coordi- 
nate q to the rate of dissipation in the system when a gen- 
eralized force acts directly on this coordinate. More specif- 
ically, the thermal noise spectrum of q at temperature T is 
given by [13] 



Saif) 



Re[Z(/)] 



(43) 



where / is frequency, Z(/) is the mechanical impedance (in- 
verse of admittance), or 



Z(/) = -2mfq{f)/F{f) 
Alternatively, imagining a sinusoidal force 

F{t)=FQCos,{2nft) 



(44) 



(45) 



with amplitude Fq acting directly on q, Eq. (43) can also be 
written as 



4k bT IVdiss 4^Br u 



F2 

^0 



p2- 
^0 



(46) 



where Wdiss is the energy dissipated per cycle of oscillation 
divided by 2n (in other words, Wjiss is the average energy loss 
per radian), U is the peak of the stored energy in the system, 
and (p is the loss angle, defined by 



(/)=Re[Z(/)]/Im[Z(/)] 



(47) 



It is important to note that (j) is in general frequency depen- 
dent. However, for an elastic body, if the frequency is low 
enough (much below the first eigenfrequency), then U can be 
computed using the quasi-static approximation, because it is 
equal to the elastic energy stored in the equilibrium configu- 
ration when a constant Fo is applied to the system. 
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B. Mechanical Energy Dissipations in Elastic Media 

It is straightforward to apply Eq. (46) to calculate the ther- 
mal noise component due to fluctuation of the position of the 
coating-air interface — the weighted average [c.f Eq. (35)] of 
the first two terms of Eq. (13). This can be obtained by apply- 
ing a force F with a pressure profile proportional to /(.?) on to 
the mirror surface (coating-air interface). In this case, elastic 
energy can be divided into bulk energy Ub and shear energy 
Us [Chapter I of Ref. [14]], with 



coating 



Ub + Us 



coating 



K 



-02+^E,-,.E,-,h/y, (48) 



where is the expansion, and E,,- is the shear tensor (see Ap- 
pendix B for details). If we give small imaginary parts to K 
and H, writing 

K^K{l+ic^B), il = fi{l+i(^s) (49) 
then Wdiss can be written as 

Wd,ss^^BUB + (^sUs (50) 

Here have introduced the loss angles (j)B and (j)s, which are 
associated with the dissipation of expansion energy density 
and the shear energy density, respectively. Note that our way 
of characterizing loss differs from previous work by Harry, 
et. al. [11], because for isotropic materials, 0b and 0^ are the 
two fundamentally independent loss angles that characterize 
the dissipation of bulk and shear elastic energy; were we to 
literally adopt and 0|| as done in Ref. [11], there would 
be modes of external driving that lead to negative dissipative 
energy, as shown explicitly in Appendix C. 

Once we have introduced (pQ and (ps, other elastic moduli 
also gain small imaginary parts correspondingly. For exam- 
ple, for the most widely used Young's modulus and Poisson 
ratio, because 



K 



3(1-2(7)^ 



2(1- 



we can write 



with 



and 



? = F(l+/0y) 

(1 - 2a)^B + 2(1 + a)^s 



a = CJ+-(l-2(j)(l+(j)(0B-0j). 



(51) 



(52) 



(53) 



(54) 



Since - 1 < (J < 1 /2, we have ( 1 - 2(j) ( 1 + a) > 0, therefore 
a has a positive imaginary part as 0b is greater than (ps, and 
vice versa. To understand the physical meaning of the imagi- 
nary part of Poisson ratio, one has to realize that Young's mod- 
ulus and Poisson ratio together describe the elastic response 
of a rod. Suppose we apply an oscillatory tension uniformly 



along a rod at a very low frequency, whether the area of the 
rod leads or lags the length of the rod depends on the relative 
magnitudes of the bulk and shear loss angles. In the situation 
when the two loss angles 0b and (ps are equal to each other, the 
Poisson's ratio is real, and we only need to deal with one loss 
angle 0^ — although there is reason to assume the equality of 
these two angles. 

If the coating material is made into the shape of a one- 
dimensional rod, and if we only consider its elongational, 
bending or torsional modes, then the Young's modulus is the 
appropriate elastic modulus associated with these modes, and 
0^ is the appropriate loss angle to apply. However, this is not 
directly relevant for coating thermal noise. An elastic modulus 
that will actually prove useful is that of the two-dimensional 
(2-D) flexural rigidity of a thin plate made from the coating 
material. 



D 



D|(1+/0d) 



12(1 -ct2; 

where h is the thickness of the plate, with 



0D = 



(1 - g - 2g2)0B + 2(1 - g + a^)<j)s 
3{l-a) 



(55) 



(56) 



As we shall see in Sec. VII A, this D is most easily measured 
through the quality factor of drum modes of a thinly coated 
sample — although this will not turn out to be the combination 
of loss angle that appear in the thermal noise. 



C. Thermal Noise of a Mirror Coated with one Thin Layer 

In the case where the coating thickness is much less than the 
size of the mirror substrate and the beam spot size, the elastic 
deformation of the substrate is not affected by the presence of 
the coating. As a consequence, if we include the elastic energy 
stored in the substrate f/sub with loss angle 0sub, we can write 



Wdiss — 0subt^sub^ 
~ 0sub + 0i 



0B% 



<psUs 
Us 



sub 



05 



sub 



£^sub 



(57) 



With the assumption of thin coating and half-infinite substrate, 
the total strain energy stored in the sample can be considered 
as t/sub- In such a way the coating adds on to substrate loss 
angle as additional, effective angles 



Ub Us 

0coated = 0sub + JJ—'pB + TJ—<l>S 
t^sub t^sub 



(58) 



Note that when the total coating thickness I is much less than 
beam spot size wq, we have t/s/t/sub ^ Us/U^uh ^ I/wq 1. 
Unfortunately, however, 0b and 05 are found to be so much 
larger than the substrate loss angle 0sub that in practice coating 
thermal noise still dominates over substrate thermal noise. 

Now suppose we would like to measure a weighted average 
of the position of the mirror surface. 



1' 



d^xw{x)z(x) 



(59) 
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with [Cf. Eq. (37)] 



w(x) ■ 



(60) 



and z{x) the position of the coating-air interface at transverse 
location x. 

According to Sec. Ill A, we need to apply a pressure profile 

of 



(61) 



onto the upper surface of the coating, which we shall also refer 
to as the coating-air interface. Straightforward calculations 



give 



Ub 

^0 



(1-2(7,)/ 



Y, (l-2g,)^(l + (7,)^ 

Y? (1-CT,)2 

1 2(l-2g,)(l + g,)(l + g,) 



1 (1 



\2-\ 



J w^{x)d^x 



^0 



21 

y 



Y, (l-g, + g,2)(l+g,)2(l-2g,)2 
72 (l-ae)2(l + a,) 

(l + (j,)(l-2g,)(l-2g,)(l + g,) 

Ys{l~cjc)^ 
(l-2g,)2(l + g,) 
y,(l-a,)2 



(62) 



J w^{x)cf'x 



(63) 



Here I is coating thickness; for Young's modulus Y and Pois- 
son's ratio (7, substrates c and s represent coating and sub- 
strate, respectively. Directly following Eqs. (46) and (50) will 
give rise to a noise spectrum of 



4^ 



Us 





(64) 



where Ub/Fq and Us/Fq are given by Eqs. (62) and (63) re- 
spectively. 

Here we can define 



J w^{x)d^x- 



Jd^xfjx) 
[Jd^xl{x)f 



1 



(65) 



as the inverse of an effective beam area. Therefore noise 
power in q is proportional to coating thickness and inversely 
proportional to beam area. In particular, for a Gaussian beam 
with 



/ 2x2 
I{x) exp ^ 



(66) 



the effective area is jz^ff = ttvvq. 

Let us compare our results to previous calculations using 
and 0||. As it turns out, if we assume (ps — 0^, then for- 
mulas for thermal noise agree with Eq. (22) in Ref. [11]. To 



illustrate the different roles now played by (pg and (ps, let us 
take the very simple case of F = = Ys and a — Gc = d?, 
where 



8Ub 41 



dUs 21 



F(f 3F^eff 



(1+(J)2(1-2(J) (67) 



F(2 3F£/eff 



(l+cj)(l-2(j)2 



(68) 



Using Eq. (64), we can get the power spectral density of the 
single layer non-penetration coating thermal noise as 



8/tBr(l-cJ-2(j2)/ 
^nfY^ 



[2(l+a)fe + (l-2(j)05]. (69) 



From Eq. (69), we can see that the bulk loss and shear loss 
contribute differently to the total noise. More importantly, at 
least in this very simple case of = Ys, the combination of (pB 
and (ps, approximately 20^ + (ps, that enters the thermal noise 
apparently differs significantly from the combination 0tot ~ 
(pB + 2(j)s, which has been measured by ring-down experiments 
that have been performed so far [15-17] — as we will see in 
Eq. (110) and will be discussed in detail in the rest of Sec. VII. 



D. Discussions on the correlation structure of thermal noise 

Before proceeding to more detailed calculations of Brown- 
ian noise that involve light penetrating into the coating layers, 
we would like to gain more insight about thermal noise by 
inspecting our existing expressions of coating thermal noise 
[Eqs. (62)-(64)] more carefully. We note that 



J w^{x)d^x. 



(70) 



where the coefficient of proportionality depends only on ma- 
terial property. From such a dependence on coating and beam 
geometries, we deduce that (i) each point on the coating-air 
interface fluctuates along the z direction independently, and 
(ii) materials at different z's within the coating also contribute 
independently to coating thermal noise. These observations 
will be confirmed mathematically in the next section. 

Finally, within the coefficient of proportionality 
[Cf. Eqs. (62) and (63)], we found three types of de- 
pendence on the Young's moduli of the coating and substrate 
materias: terms proportional to l/F, are expected to arise 
from fluctuations in coating thickness, terms proportional 
to Yc/Y} can be interpreted as arising from coating thermal 
stresses driving the substrate-coating interface, while terms 
proportional to 1 /Ys are therefore interpreted as correlations 
between the above two types of noise. 



IV. CROSS SPECTRA OF THERMAL NOISE 
COMPONENTS 

In this section, we compute the cross spectra of each com- 
ponent of coating thermal noise, and assemble the formula 
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for the spectral density of the total noise. Specifically, in 
Sec. IV A, we compute the cross spectra of the thickness fluc- 
tuations any two uniform sublayers of the coating, and obtain 
the cross spectrum of S^^; in Sec. IV B, we compute the cross 
spectra involving height fluctuation Zs of the coating-substrate 
interface, i.e., Ss-z^ and 57^,^; in Sec. IV C, we dissect the 
above results and analyze the separate roles of bulk and shear 
fluctuations; in Sec. IV D, we write down the full formula of 
coating thermal noise. 



A. Coating-Thickness Fluctuations 

Let us start by calculating thickness fluctuations of individ- 
ual layers and correlations among them. Following Levin's 
approach, we imagine applying two pairs of opposite pres- 
sure. 



/i (x) = -Pb w 1 (x) , /3 (x) = Fqwj, (x ) 



(71) 



in the z direction on layer I and layer III, as shown in Fig. 2, 
with thickness of li and Z3, respectively. Here wi{x) and 
W3(x), like the w{x) used in Eq. (59), provides the shape of 
the pressure profiles. 

We assume that within each of I and III, there is only one 
type of material, yet there could be arbitrary number of differ- 
ent material sub layers in II. As it will turn out, the precise lo- 
cations of layers I and III along the z direction does not affect 
the result, as long as they do not overlap, or in other words, 
layer II has non-zero (positive) thickness. 

Throughout this paper, we shall assume that the beam spot 
size is much less than the radius of the mirror, so that we 
can make the approximation that the mirror surface is an infi- 
nite two-dimensional plane. In this case, we perform a spatial 
Fourier transformation for the applied pressure. 



(72) 



f,{k) = / e''-'fj{x)d'x = FQWj{k) , - 1,3, 



and carry out our calculations for strain and stress distribu- 
tions in the coating-substrate system in the Fourier domain. 

We further assume that the coating thickness is much less 
than the beam spot size, which is inverse the maximum spa- 
tial frequency contained in m>i 3. This means we only need to 
consider ^'s with ^ 1, with I the total coating thickness. 
According to calculations of Appendix B, non-zero compo- 
nents of the stress and strain tensors in Layers I and III are 
found to be (in the spatial Fourier domain) 



(JlM>l 
1-C71 



fi, = w\Fq, 



-1 (l-2gi)(l+gi)vvi ^ 

= TTVl ^ rQ , 



and 



fm 

XX "yy 



fill _ ^3^-3 p .^01 _ p 
1 — C73 



(73) 



(74) 



(75) 



cIII 



(l-2g3)(l+g3)H'3 
¥3(1-03) 



(76) 



respectively. 

Note that deformations within layer I only depends on wi 
(not ^3), while deformations within layer III only depends on 
W3 (not wi) — while regions outside these layers are found 
to have vanishing strain and stress. This means we can treat 
deformations caused by each pair of forces independently, as 
long as layer I and layer III do not overlap. The deforma- 
tions are also independent from the thickness of the layers. 
The vanishing of deformations outside these layers means that 
when we introduce additional pairs of opposite forces, the new 
deformations introduced will be constrained within those new 
layers — as long as those new layers do not overlap with ex- 
isting ones. This independence originates from the linearity of 
elastic response, and the fact that coating strains induced by 
force applied on a single surface within the coating, as dis- 
cussed in Appendix. B, does not depend on distance away 
from that surface, as seen in Eqs. (B25)-(B32). The situa- 
tion here is analogous to the electric field generated by several 
pairs of oppositely-charged infinite parallel planes. 

In terms of thermal noise, such a distribution of elastic de- 
formations corresponds to a dissipation energy that consists of 
two independent terms, each corresponding to one layer and 
proportional to its thickness: 



f2 



Wnh J w]d^x + W33l3 J wjd^x (77) 



Here we have defined, for J ~ 1,3: 



W 



J J - 



_ ( l-2a,)(l + g,) 
3{l -aj^Yj 



^1 



0i + (l-2(T,)</)i 



This means the fluctuation of 



1 = j [wi{x)5li{x) + W3{x)5l3{x)\d'^x 



is given by 



y 



Wjjlj I wjix)d'x 



(78) 



(79) 



(80) 



The absence of cross terms here means that fluctuations in 
5li{x) and 5l3{x') are uncorrected — and hence statistically 
independent. Furthermore, within each layer, in the same 
spirit as the discussions in Sec. HID, the particular form of 
dependence on Ij and Wj{x) indicates that Szz fluctuations at 
different 3-D locations (within this layer) are all uncorrected 
and have the same spectrum. In this way, we obtain the cross 
spectral density of Szz two arbitrary 3-D locations within 
the coating: 



{x,z;x',z') 



4kBT 



5,j8^^\x-^)5{z-z')Wjj (81) 



Here we have assumed that {x,z) belongs to layer /, while 
{x,z) belongs to layer j. (The association to layers helps to 
identify the material property to be used in Wjj.) 
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Figure 2: Illustrations of forces applied onto various interfaces within 
the coating. Each of Layers I and III in the coating are assumed to 
be uniform (but they might each contain a different kind of mate- 
rial); region II denotes the entire gap between them, which may well 
contain many different dielectric layers. A pair of force distribution 
/l (/s) in opposite directions is exerted on opposite sides of Layer I 
(III), while fs is exerted on the coating-substrate interface. The three 
distributions may well have different profiles (as also illustrated in 
the figure). 



B. Fluctuations of Coating-Substrate Interface and their 
correlations with coating thickness 

To investigate the correlation between height of the coating- 
substrate thickness, Zs{x) and the thickness of each coating 
layer, 5lj{x), we apply a pair of pressures f\ (x) = Fqwi (x) at 
opposite sides of Layer I, and force fs{x,y) — FqWsQc) onto the 
coating-substrate interface (along the —z direction), as shown 
in Fig. 1. The same strain and stress as in Eqs. (73) and (74) 
are driven by fi , which are only non-vanishing within layer I. 
On the other hand, drives uniform strain and stress over the 
entire coating, with non-vanishing components of stress and 
strain given by. 



|7;-,II = 



-2(7,2)7, 



(l + a,)fc2y, 



kxky 




1— (jj. 





(82) 



where Young's modulus and Poisson's ratio a, of the coat- 
ing are given by values within layer I. The total dissipation in 
this case will have the following structure. 



p2 
^0 



Wn J w\cf-x + 2Wis j 



w\Wsdx- 



Wss I wl(fx 



(84) 

with the first term arising from dissipation in layer I that is due 
to strain and stress driven by f\ , the second term also arising 
from dissipation in layer I arising from cross terms between 
strains and stresses caused by f\ and and the third term 
arises from dissipations throughout the entire coating, due to 
strain and stress caused by . Here Wn is the same as defined 
by Eq. (78), and 



;i-CJ,-2(j2)(i_(j,-2(j2) 



2(l-a/)2y, 



(l-(T,)2y2 



i(j>^-<j>^) (85a) 
1 



1 



(85b) 



Note that we have added a superscript ( j) for Wss to indicate 
that here the dissipation is due to force applied on one thin 
layer alone. 

Here again, the dependences on w\ and vv, indicates that 
fluctuations at different transverse locations, x ^ x', are un- 
coiTelated, while the li in front of Wn, and the arbitrariness 
of h means that S^z fluctuations at different z locations within 
the thin layers are uncoiTelated. The li in front of both Wis 
and W„ indicates that all S-,z within layer I are correlated with 
Zs the same way, even though all of them are mutually uncor- 
related. 

This allows us to extract the following 



Ss„s-{x;x,z! 



3nf 

— 5{x~x)W,. 



(86a) 
(86b) 



Here for Eq. (86b), j is the layer with which z! is associated; 
and this labeling is to help identify which material parameter 
to use in Wjs- 



C. The natomy of coating thermal noise 



w,(l-(Jv-2g2; 



^2 

kxkx' 



kxk\ 



Here let us assemble Eqs. (81), (86a) and (86b) from the 
(83) previous sections, and write: 
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,z') 



4^Br(l+(7y)(l-2(J;)rl+CJ, 



3nf Yjil-ajY 

AksT (1-CJ,- 2(7,2)2 



3nf 



YI 

.. (1-cJ,- 



1-2(7; 



1 - (7y + (7f 



1+C7,- 



2^Br(l-a,-2a/)(l-(J,-2(j/ ^ 

w ^(T^^;? '^'•'"^'^■^^ 



(87a) 
(87b) 
(87c) 



Here we have assumed that z belongs to the i-th layer and 
that z' belongs to the y-th layer, respectively. The thickness 
fluctuation of different layers are mutually independent [note 
the Kronecker delta in Eq. (87a)], while thickness fluctuation 
of each layer is correlated with the height fluctuation of the 
coating-substrate interface [Eq. (87c)]. 

Fluctuations described by Eqs. (87a)-(87b) can be seen as 
driven by a set of microscopic fluctuations throughout the 
coating. Suppose we have 3A^ thermal noise fields (i.e., 3 for 
each coating layer), «y (x), n^j'^ (x) and n^^ (x), all independent 
from each other, with 



(88a) 



37tfYj{l-Gj)^ 



(88b) 



and all other cross spectra vanishing. Here /' labels coating 
layer, the superscript B indicates bulk fluctuation, while Sa 
and Sb label two types of shear fluctuations. The normaliza- 
tion of these fields are chosen such that each of these fields, 
when integrated over a length Ij along z, have a noise spectrum 
that is roughly the same magnitude as a single-layer thermal 
noise. 

Noise fields (x), ny'*(x) and can be used to generate 
thickness fluctuations of the layers and the interface fluctua- 
tion (87a)-(87b) if we define 



u,,{x,z) = Cy'j{x,z)+Cfnf (x,z) 



(89) 



and 



f J 



D'Jn']{x,z)+Dfnf{x,z 



Dfnf{x,z) 



(90) 



For each coating layer, and Z)^ are transfer functions from 
the bulk noise field to its own thickness 51 j and to sur- 
face height Zs, respectively; and are transfer func- 
tions from the first type of shear noise to thickness and sur- 



face height; finally is the transfer function from the sec- 
ond type of shear noise to surface height (note that this noise 
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Figure 3: Illustration of the correlations between coating thickness 
81 j and the height of the coating-substrate interface, Zs- On the left, 
for a bulk deformation: when a coating element is expanding, its 
expansion along the x-y plane lifts the coating-substrate interface up- 
wards, causing additional motion of the coating-air interface corre- 
lated to that caused by the increase in coating thickness. On the right, 
a a particular shear mode: when a coating element is expanding, its 
contraction along the x-y plan pushes the coating-substrate interface 
downwards, causing addition motion of the coating-air interface anti- 
correlated to that caused by the increase in coating thickenss. 



Thickness (Sj) 



Surface height izs) 



Bulk : 

Shear- Sa 

A j 

Shear 
B 



D 



I -a,- laj Yj 
1 - CT, - 20- 2 Yj 



(none) 



2^- 
^/3(l-CT,■)(l- 



2oj Ys 



2^1-2(7^-(l + (Tj) Y, 



Table I: Transfer functions from bulk and shear noise fields to layer 
thickness and surface height. 



field does not affect layer thickness). Explicit forms of these 
transfer functions are listed in Table. 1. 

Equations (89) and (90) owe their simple forms to the un- 
derlying physics of thermal fluctuations: 

For bulk noise, i.e., terms involving n^, the form of 
Eqs. (89) and (90) indicates that the interface fluctuation due 
to bulk dissipation is simply a sum of pieces that are directly 
proportional to the bulk-induced thickness fluctuations of each 
layer This means the thermal bulk stress in a layer drive si- 
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multaneously the thickness fluctuation of that layer and a fluc- 
tuation of the coating-substrate interface. The fact that and 
Cj having the same sign means that when thickness increases, 
the interface also rises (with intuitive explanation shown in 
Figure 3). This sign of correlation is generally unfavorable 
because the two noises add constructively towards the rise of 
the coating-air interface. 

For shear noise, the situation is a little more complicated, 
because unlike bulk deformations, there are a total of 5 pos- 
sible shear modes. From Eq. (73) and (74), it is clear that 
/i , applied on opposites of Layer 1 (Figure 2), only drives the 
xx+yy — 2zz shear mode and the xx+yy+ zz bulk mode, while 
from Eq. (82) and (83), the force distribution drives three 
shear modes of xx — yy, xy + yx, and xx+yy ~ 2zz- This means 
while thermal shear stresses in the xx + yy — 2zz mode drives 
layer thickness and interface fluctuation simultaneously, there 
are additional modes of shear stress, xx — yy and xy + yx, 
that only drives the interface without driving layer thickness. 
Our mode Sa, which drives both layer thickness and inter- 
face height, therefore corresponds to the physical shear mode 
of xx + yy — 2zz', our mode Sb, which only drives interface 
height, corresponds to the joint effect of the physical shear 
modes xx — yy and xy+yx. It is interesting to note that for Sa, 
its contributions to 81 j and Zs are anti correlated, because C'^" 
and D'^i have opposite signs. This is intuitively explained in 
Fig. 3. 

As an example application of Eqs. (89) and (90), if we ig- 
nore light penetration into the coating layers, namely, when 
thermal noise is equal to 



;np - 



(91) 



D. Full formula for thermal noise 

As we consider light penetration into the coating, we resort 
to Eq. (24), and write: 



^{x)-iax) 



L dz 



iSjiz) 



+ 



2 



nj{x,z) 
nf{x,z) 



(93) 



Here spectra of independent fields n^, n-^ and have been 
given in Eqs. (88a)-(88b), e is defined in Eq. (25), while the 
transfer functions C's and D's are listed in Table. 1. 

We can then obtain the spectrum of phase noise (after aver- 
aging over the mirror surface, weighted by the power profile 
of the optical mode) as 



pi dz 
f'j dz 



1 -Im 



1 -Im 



2 



S] 



D 



2 Ij ,S 



(94) 



we have 



:np 



dz 



-{d/+Df)nf 



(92) 



in which contributions from each layer has been divided into 
three mutually uncorrected groups, each arising from a dif- 
ferent type of fluctuations. Here we see explicitly that and 
sharing the same sign increases contributions from rfi, C^^ 
and D^'^ having opposite signs suppresses contributions from 

Finally, we note that in the spectral density of ^"P, contri- 
butions directly from coating thickness will be proportional to 
|Cj*p and |C^'*p, and hence proportional to l/Yc, those from 
interface height will be {ofl"^ and \D^f\2, and hence 

proportional to Yc/Y^, while those from correlations will be 
proportional to C^D^ and CJ'*D^'', and hence proportional to 
1 /Fj. This confirms our anticipation at the end of Sec. Ill D. 



and spectrum of amplitude noise as 



f^j dz 



n 2 



(95) 



Here A; is the wavelength of light in Layer j, and we have 
defined 



5J: 



4kBTXj^l{l-aj-2aj) 
3nfYj{l-ajy£/,s 



X=B,S. 



(96) 



which is at the level of coating thickness fluctuation of a sin- 
gle layer of dielectrics with material parameters identical to 
layer j and length equal to Xj. Note that the quantity S^ only 
depends on the material properties (and temperature) of the 
layer, and is independent from length of that layer; the quan- 
tities qj, on the other hand, give us the relative thermal-noise 
contribution of each layer in a dimensionless way. 

Note that the reason for keeping the integrals in Eqs. (94) 
and (95) is because £ has a z dependence, which originates 
from the fact that the back-scattering contributions to 5(p/s 
and dr/sa weighted integral of u„ within each layer [Cf. (17) 
and (18)]. 
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Figure 4: Two basic transformations involved in solving for optical 
fields in a multi-layer coating. 



Parameter 



Tantala(Ti205) Silica(Si02) 



Refractive index 
Poisson's ratio 
Young's modulus (Pa) 
Loss angle (fe = 0^) 
Photoelastic coefficient 



2.07 [18] 
0.23 [19] 
1.4 X 10" [20] 
2.3 X 10"'*[21] 
-0.50 [23] 



1.45 [18] 
0.17[19] 
7 X 10i°[19] 
4.0 X 10"^ [22] 
-0.41 [24] 



Table II: Baseline material parameters. 
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e-'<l> 




b 



(99) 



In this way, assuming the input and output field amplitude 
at the top surface of a multi-layer coating to be v'l and V2, and 
writing those right inside the substrate to be si and S2, we have 



(100) 



(101) 



(102) 
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V2 




V2 



where M is given by 

The complex reflectivity is given by 

p = -M21/M22 



V. EFFECT OF LIGHT PENETRATION INTO THE 
COATING 

In this section, we synthesize results from Sec. 11 and 
Sec. IV, and compute the full Brownian thermal noise for 
coating configurations. We will illustrate how the light pene- 
tration affects the total noise in highly reflective coatings. 



A. Optics of multi-layer coatings 

For completeness of the paper, we briefly review how light 
penetration coefficient d log p/d(j)j can be calculated. 

From an interface from layer / to j (here j is either / + 1 or 
we denote the reflectivity and transmissivity of different 



layers by r,y and tij: rf 



ri; + n i 



(97) 



We also define hn+i =n\, since that is the refractive index of 
the substrate. 

A matrix approach can be applied to solve for the amplitude 
of light inside the layers, when we view the coating as made 
up from two elementary transformations, each representable 
by a matrix. In this approach, instead of writing out-going 
fields in terms of in-going fields, one write fields to the right of 
an optical element in terms of those to the left. As illustrated 
in Figure 4, for reflection at an interface (left panel), we write 



(98) 



On the other hand, for propagation across a gap with phase 
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B. Levels of light penetration in Advanced LIGO ETM 
Coatings 

In Advanced LIGO, the coating stack is made from two 
alternating layers of materials: Si02 (ni = 1.45) and Ta205 
(n2 — 2.07). Here we consider the End Test-mass Mirror 
(ETM). In order to achieve very high reflectivity, the coating 
is made of 19 successive pairs of alternating Si02 and Ta205 
layers, all A/4 in thickness except the top one, which is A /2. 
We will refer to this as the conventional coating. An alterna- 
tive design has been made to allow the coating to operate at 
both 1064 nm and 532 nm. We shall refer to this as the Ad- 
vanced LIGO coating (see Appendix. D) [25]. 

In Fig. 5, we plot real and imaginary parts of dlogp /d^j 
and dlogp /drj, for both conventional and Advanced LIGO 
coating. Here we note that the real parts of these derivatives 
are at the order of 10^^, which means ^ is less than | by 
6 orders of magnitude. This, together with considerations in 
Sec. II E, will make amplitude coating noise negligible. 

In Eq. (27), we have divided contributions to ^ into four 
terms, the first, Zs, is the height of the coating-substrate inter- 
face, while the other three are related to fluctuations in layer 
thickness, 8lj, 51'j and 5lj, see Eqs. (27)-(30). We can il- 
lustrate the effect of light penetration by showing the relative 
size of these three contributions for each layer. In Figure 6, we 
carry out this illustration, for conventional coating on the left 
panel and for Advanced LIGO coating on the right. We use 
solid black line to indicate the non-photoelastic part of 
[i.e., term not containing jSy, see Eq. (28)], and we use red- 
long-dashed, blue-short-dashed, and purple-dotted curves to 



indicate the photoelastic part of /{{dljf) 
{{Sl^jY) / {{5ljY), respectively. The weighting 



and Sr^' 
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Figure 5: Real (solid curves) and imaginary (dashed curves) parts of d logp/d(j)j (upper panel) and dlogp/drj (lower panel), for conventional 
(red curve) and Advanced LIGO (blue curve) coatings. [Note that Re{d log p/d(j>j) = for conventional coating.] 



factors. 



C. Thermal noise contributions from different layers 



(103) 
(104) 



have been added for 5^''' and •S^f'', respectively, to correct for 
the fact that Slj and 5lj have different r.m.s. values compared 
to 51. Because of the lack of experimental data, we have as- 
sumed j5j — —0.4 identically. Note that in order to focus on 
the effect of light penetration, we have only showed the first 
10 layers. 

In the figure, the effect of light penetration into the coating 
layers is embodied in the deviation of the black solid curve 
from unity in the first few layers, and in the existence of the 
other curves. Although we cannot perceive the correlation be- 
tween these contributions, we can clearly appreciate that only 
the first few layers are penetrated, and that the total effect of 
light penetration will be small. We should also expect the ef- 
fect of photoelasticity (dashed curves) to be small, and the 
effect of back-scattering (which gives rise to ,3^.^' and \ 
blue and purple dashed curves) even smaller. 



Let us now examine the breakdown of the total coating 
noise by plotting the coefficients and q^j in Eq. (94). In 
Fig. 7, we plot silica contributions on top panels, and tan- 
tala contributions on lower panels, with bulk contributions on 
left panels, and shear contributions on right panels. Here we 
use the baseline parameters shown in Table 11. As it turns out, 
the results for conventional and Advanced LIGO coatings are 
hardly distinguishable from each other — therefore we only 
use the Advanced LIGO coating. The red curve uses )3 = — 1, 
black uses j3 = and blue uses P — I- Superimposed onto the 
solid lines are dashed lines of each type, calculated without 
introducing the back-scattering terms; the effect is noticeable 
for the first few layers. 



VI. DEPENDENCE OF THERMAL NOISE ON MATERIAL 
PARAMETERS 

Experimental knowledge of coating materials is limited. 
Most notably, values of Young's moduli and Poisson's ratios 
of the coating materials are still uncertain, while only one 
combination of the two loss angles have been experimentally 
measured by ring-down experiments. In this section, we ex- 
plore the possible variation in coating Brownian noise, away 
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Figure 6: Light penetration into the first 10 layers of a 38-layer coating (left panel for conventional coating and right panel for Advanced 
LIGO coating). We plot the non-photoelastic part of in black sold curves, the photoelastic part of in long-dashed red curves, as well as 
.3^*(scaled by rms value of Sl'j with respect to the rms value of Slj, shown in short-dashed blue curves) and T? (scaled by rms value of Slj, 
shown in dotted purple curves). These plots indicate that for both structures, light penetration is restricted within the first 10 layers. 



from the baseline configuration (Table 11), considering these 
uncertainties. We shall use the Advanced LIGO coating struc- 
ture mentioned in the previous section. 



A. Dependence on Ratios Between Loss Angles 

In the baseline (Table II), we have assumed that (j)B and 
05 are equal, but this is only out of our ignorance: experi- 
ments have only been able to determine one particular combi- 
nation of these two angles. We now explore the consequence 
of having these loss angles not equal, while keeping fixed the 
combination measured by ring down rate of drum modes [see 
Eq. (??)]. 

In Figure 8, while fixing all other baseline parameters, we 
plot how each type of thermal noise (i.e., silica vs tantala, bulk 
vs shear) varies when the ratio (j)B / (j)s for both tantala and sil- 
ica layers varies between 1/5 and 5. We use blue for tantala, 
red for silica, dotted for bulk, dashed for shear, and solid for 
the total of bulk and shear In this configuration, tantala lay- 
ers' contribution to thermal noise always dominate over silica 
layers, mainly due to the higher loss angle. As we vary the 
ratio between the loss angles, there is moderate variation of 
thermal noise. For the dominant tantala, as ^b/^s vary from 
1/5 to 5, there is a 30% change in thermal noise, while for 
silica, the change is a more significant 68%. 

As we see from Fig. 8, a larger value of 0b /^s gives rise to 
higher bulk, lower shear, and higher total noise — this is rea- 
sonable because bulk fluctuations drive correlated noise be- 



tween layer's thickness and the height of coating-substrate in- 
terface, while shear fluctuations drive anti-correlated noise, as 
shown in Fig. 3. 

Moreover, the fact that variation is more significant for sil- 
ica layers can be explained when we recall that thickness- 
induced thermal noise is proportional to \/Yc, while surface- 
height-induced thermal noise is proportional to Y^-jY'^. For sil- 
ica layers, Yc is assumed to be equal to F^, so the two types of 
noise being added (bulk) or subtracted (shear) are more com- 
parable in magnitude; by contrast, the Young's modulus of 
tantala layers is significantly higher than that of the substrate, 
causing the noise to be dominated by fluctuations of the height 
of the coating-substrate interface, therefore making correla- 
tions between the two types of noise less important. 

In Fig. 9, we plot variations in the total noise as we vary 
/ 05 for silica layers (blue) or tantala layers (red) only, and 
fix the other one. It shows that the variance of tantala's loss 
angle will generate larger change of the total noise. 



B. Dependence on Young's moduli and Poisson's ratios 

Since the Young's modulus and Poisson's ratios of coating 
materials, especially of tantala, are also uncertain. In Fig. 10, 
we plot variations of tantala thermal noise when its Young's 
modulus varies from the baseline value by up to a factor of 2, 
for 0b/05 = 0.2, 0.5, 1, 2 and 5. The noise is seen to vary by 
^15% as Young's modulus varies by a factor of ^ 2. 

We can also explain the way the thermal noise varies as a 
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Figure 7: A break-down of thermal noise contributions from silica (upper panels) and tantala (lower panels) layers, from bulk (left panels) and 
shear (right panels) losses. Blue curves correspond to j3 = — 1, black j3 = and red jS = 1. Dashed curves indicate results calculated without 
including back-scattering effects. 



function of F^.. Starting from the baseline value, a lower 
leads to a lower thermal noise, until Yc becomes comparable 
to Ys (which we fix at the baseline value, equal to O.Sira), and 
starts to increase again. Such a behavior is reasonable be- 
cause thickness noise spectrum and interface noise spectrum 
are proportional to ^ 1 /Yc and ^ Y^/Y^, respectively — as 
we decrease Yc from the baseline ixa value, we transition from 
interface fluctuation being dominant towards equal amount of 
both noises (which gives a minimum total noise), and then 
towards thickness fluctuation becoming dominant. 

In Fig. 1 1, we explore the effect of varying coating Pois- 
son's ratio, for the same values of ^b/^s chosen in Fig. 10. 
In the baseline assumption of 0b = 0s, when bulk and shear 
have the same level of loss, thermal noise does not depend 
much on Poisson's ratio. However, if 0b/0s turns out to differ 
significantly from 1, and if Poisson's ratio can be larger than 
the baseline value by more than ~ 0.1, then thermal noise can 
vary by ^ 10%. 



vary j3 between -1 and + \. Dashed curves are obtained ignor- 
ing back-scattering effects. 

It is interesting to note that for small values of j3, the depen- 
dence of noise on /3 have different trends for bulk and shear 
contributions. This is also related to the different types of cor- 
relations between thickness and interface height fluctuations. 
As we can see from the Figure, the effect of varying )3 is small, 
since it only affects thermal noise due to light penetration into 
the first few layers. If bulk and shear losses are indeed compa- 
rable, then cancelation between these two types of noises (es- 
pecially for the more lossy tantala layers) will likely make the 
photo elastic effect completely negligible. Even in the case 
when one particular type of loss dominates shall we expect 
at most ^2% contribution from photo elasticity of the more 
lossy tantala — if we further assume that |j3 1 ~ 1 [right panel 
of Fig. 12]. 



D. Optimization of Coating Structure 



C. Dependence on Photoelastic Coefficients 

Photoelastic properties of the coating materials are not yet 
well known. In Fig. 12, we plot the fractional change in ther- 
mal noise, separately for silica (left panel) and tantala (right 
pane), and for bulk (blue) and shear (red) losses, when we 



Although a standard highly reflective coating consists of 
A /4 layers of alternating material capped by a A/2 layer, this 
structure can be modified to lower thermal noise while still 
maintaining a high reflectivity for the 1064nm carrier light, 
e.g., as shown by Agresti et al. [26]. As their results have in- 
dicated, for baseline coating parameters and neglecting light 
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Figure 8: (Color Online) Variations in thermal noise contributions 
when (Pb/'Ps is varied. Contributions from tantala layers is shown 
in blue, those from silica layers are shown in red. The total thermal 
noise is in black. Bulk contributions are shown in dotted curves, 
while shear contributions are shown in dashed curves. 
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<t>B/<t>S 

Figure 9: (Color Online) Variations in total noise when ^b/^s is 
varied: (solid) total noise, (dotted) total bulk noise, (dashed) total 
shear noise. The red (blue) curve corresponds to only varying (ps/'I's 
for tantala (silica). With (ps/'Ps of tantala or silica varying from 0.2 
to 5, the change in total noise is 58.1% and 10.6% respectively. 



penetration into the coating layers [11], the optimal structure 
is more close to a stack of pairs of A/8 (Ta205) and 3A/8 
(Si02) layers, capped by a A/2 (Si02) layer. This alternative 
coating structure shortens the total thickness of the more lossy 
tantala layers, while maintaining a high reflectivity for the 
light. The Advanced LIGO type coating given in Appendix D, 
on the other hand, has been optimized considering reflectiv- 
ity at both 1064 nm and 532 nm, as well as thermal noise — 
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Figure 10: Thermal noise contribution from tantala, as its Young's 
modulus deviates from baseline value, for (j)B/(j)s=5 (blue dashed), 2 
(blue dotted), 1 (black solid), 1/2 (red dotted), and 1/5 (red dashed). 
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Figure 11: Thermal noise contribution from tantala, as its Poisson's 
ratio deviates from baseline value, for (j)B/(j)s=5 (blue dashed), 2 
(blue dotted), 1 (black solid), 1/2 (red dotted), and 1/5 (red dashed). 

although light penetration into the layers have not been con- 
sidered. 

In this section, we carry out a numerical optimization taking 
penetration into account. We first fix the number N of layers 
(N is even, so we have N/2 pairs), and then for N, we use the 
Lagrange multiplier method to search for the constrained min- 
imum of Sth, fixing Tio(4 and 7532, namely the power transmis- 
sivity, 1 — IpP assuming the coating is lossless, evaluated at 
1064nm and 532 nm, respectively. The quantity we seek to 
minimize (or, the cost function) is 

y =y^ + All rio64 + Al2 (7532 -5%)2 (105) 

As we vary fii and IJ.2 and minimizing y, we obtain the con- 
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Figure 12: Fractional change in the contribution to thermal noise from all silica layers (left panel) and all tantala layers (right panel), due to 
bulk (blue) and shear (red) loss. Dashed lines indicate results calculated without including back-scattering terms. 



strained minimum of \/5jh for different pairs of (r532, 7\o64)- 
The aim is to obtain a series of coating configurations with ap- 
proximately 5% transitivity for 532 nm, and with minimized 
thermal noise for a variable 3 - 20ppm transmissivity for 
1064 nm. (Note that the choice of the cost function contains a 
certain level of arbitrariness.) 

Since we are going to carry out minimization for a large 
number of multipliers over a large number of degrees of free- 
dom, we have chosen to proceed gradually allowing only the 
first n pairs and last n pairs of layers to vary, while maintaining 
the same pair structure foiN/2 — n pairs in the middle (repeat- 
ing doublets). In other words, our coating structure looks like: 



free 



repeating pair 



free 



2n layers n - 2n layers 2n layers 

In this work, we found that it suffices to choose n = 2 (which 
corresponds to optimizing over 10 parameters); further in- 
creasing n does not lead to noticeable improvements. Dur- 
ing our numerical optimization, we have adopted the downhill 
simplex method [27, 28]. 

Results for baseline material parameters (Table. II) and 
=38, 40 and 42 have been shown in Figure 13. This fig- 
ure indicates that different numbers of layers should be cho- 
sen for different target rio64 - more layers are required for 
lower transmissivity (higher reflectivity). Overall, the optimal 
thermal noise varies by around ^ 10% as for T\qm from 3 to 
20ppm. In particular, for the standard Advanced LIGO re- 
quirement of 5 ppm (see first column of Table III), 42 layers 
are found to be optimal. This is 2 more pairs or 4 more lay- 
ers than the 38-layer A/4 doublet, which has the minimum 
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Figure 13: Optimized thermal noise versus transmissivity at 
1064 nm, for a coating of 38 (red), 40 (blue), and 42 (purple) lay- 
ers. 



number of layers to reach 5 ppm. The larger number of lay- 
ers here gets lower thermal noise (by 6 %) because the more 
lossy tantala layers are shortened, and the less lossy silica lay- 
ers lengthened. 

We have further optimized the structure when the ratio 
0b/05 is different from 1, while keeping fixed the effective 
loss angle measured so far — as done in Sec. VIA. For 
rio64 — 5 ppm, we have listed results of optimized coating 
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target 
<I>b/(I>s 




Resulting Coating Structure 


/cOpt 

V th 


V th 


A' 


First 4 layers 


Repeated Pair 


Last 4 layers 


1 

~ 5 


h ~ ^ 


_ s 


1/5 


42 


0.0479 0.1581 0.3430 0.1760 


0.2919 0.1897 


0.3164 0.1738 0.3178 0.1627 


5.01 


6.64 


8.81 


5.35 


1 


42 


0.1020 0.1250 0.3267 0.1917 


0.2911 0.1914 


0.3110 0.1752 0.3196 0.1609 


5.02 


6.64 


8.81 


7.05 


5 


42 


0.1118 0.0968 0.3353 0.1882 


0.2893 0.1939 


0.3135 0.1673 0.3199 0.1662 


5.02 


6.64 


8.81 


9.33 



Table III: Results of coating-structure optimization. We list optimized coating structures for Tiq(,4 = 5ppm and 7532 = 5%, for three target 
values of ^b/^s while fixing the measured effective loss angle [Eq. (56)] and other baseline material parameters [Table II]. Thickness of 
coating layers are given in units of wavelength (for 1064 nm light). For each optimized coating structure, thermal noise is calculated separately 
for all three values of and given in units of 10^^^m/\/Hz (thermal noise for the target (j>B/<j>s is given in boldface, and boldface 

numbers should be the minimum within its column); thermal noise spectra of the 38-layer X /4 stack assuming the target tj)B / >j)s are also listed 
for comparison. 



structure and thermal noise in the second and third columns 
of Table HI. The extent of variation found here is compara- 
ble to those obtained in Sec. VIA using a standard coating 
structure without optimization: the optimal coating structures 
consistently lower thermal noise by about 6%. In addition, 
as shown in Table III, the optimal coating structure is robust 
against changes in ^b/^s'- structure obtained for any one of 
the values of the ratio is already almost optimal for all other 
ratios. 



VII. MEASUREMENTS OF LOSS ANGLES 




Figure 14: Rectangular shaped thin plate (axbx c) with thin coat- 
ing (thickness d): c <^ a,b;d <^ c. The transverse vibration mode is 
considered in this case 



In this section, we study possible mechanical ringdown ex- 
periments that can be used to measure independently the bulk 
and shear loss angles, (^b and (^5 of a coating material. 

In a ringdown experiment, a sample with a high intrinsic 
Q is coated with a thin layer of the coating material in ques- 
tion. Due to the mechanical losses in the coating, the quality 
factor of the mechanical eigenmodes of the sample will be re- 
duced [29, 30]. More specifically, for the n"' ei genmode with 
resonant frequency /„, if an e-folding decay time of T„ is mea- 
sured, then the quality factor is 

Q„ = 7zf„x„ , (106) 

while correspondingly, the loss angle is given by 

</»(./«)- l/fi,M (107) 

which is equal to the amount of energy dissipated Wdiss per 
radian. 



A. Bending Modes of a Thin Rectangular Plate 

Figure 14 shows the schematic geometry of a rectangularly 
shaped sample, in which a thin coating layer with thickness d 
is deposited on a rectangular plate with dimensions ax bx c 
(c ^ a,b), and d is much less than c. If we pay attention 
only to the transverse oscillations of the plate, the amount of 
energy stored in the coating layer, in the form of bulk and 
shear energies Ub and Us, as a fraction of the entire energy U, 



can be calculated as [Reference?] 



Ub 
U 

Us 
U 



fl'}', (1-Cj2)(l-2(J,) 



cY, (1"(T.)' 
2dY, {l-a^){l-cj, + G^) 
c y, (l-(j,)2(l + cj,) 

Using Eq. (58), the the total loss angle of the sample is 

= 0sub 



(108) 
(109) 



0b(1-2(J,)+20' 



1 - CJ, 



l+CTc 



(110) 



It is not surprising that this combination of and 0s is 
proportional to [cf. Eq. (56)], the loss angle of the 2-D 
flexural rigidity of the coating material, which we defined in 
Sec. IIIB. This is because when the drum mode of a thinly 
coated plate is excited, the stress T,, remains zero within the 
coating layer, and the layer's elastic response is governed by 
the flexural rigidity, as defined in Sec. 13 of Ref. [14]. 

As it turns out, the part of coating thermal noise due to 
bending of the coating-substrate interface [S^^^^ in Eq. (87b)] 
also depends directly on because the loss mechanism in 
this case is the same as during the oscillation of a drum mode 
— one only applies a perpendicular force from below the coat- 
ing layer, while keeping T,, = within the layer. 

It proves less straightforward to connect the thickness fluc- 
tuation part of thermal noise [S',,.,,. in Eq. (87a)] to the effec- 
tive loss angle of either Y or D. Although the loss mechanism 
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Figure 15: Thin cylindrical shell with thin coating outside. The first 
torsional eigenmodes of such a shell can be used to measure the shear 
loss angle of the coating. 



A more accurate calculation may be found by using the 
Fliigge shell theory [33]. 

Using the values from Table IV, we can estimate the res- 
onant frequency to be 9.2kHz. The coating contribution to 
loss angle, assuming a 0^ of at least 10^^, would be at least 
the order of 10^^, which seems plausible to be extracted from 
ring-down measurements. 

Table IV: Example parameters of a thin, uniformly coated cylindrical 
shell (Si02) 

L R c d 
unit(mm) 200 50 1 0.04 



here is due to the compressing of a thin membrane from both 
sides — this membrane is not characterized by vanishing T^^ 
and Tyy, because the coating is attached to a substrate which 
provides restoring forces along the transverse (x and y) direc- 
tions. However, in the case when the Poisson ratio Oc of the 
coating vanishes, the thickness fluctuation does depend on the 
loss angle of the Young's modulus. 

For our baseline parameters, mechanical dissipation is 
mostly contributed by the tantala layers, and because the 
Young's modulus of the tantala coating material is assumed to 
be much greater than that of the substrate, the largest contri- 
bution to the LIGO mirrors' Brownian noise is bending noise 
5,.,^,. This explains why the noise only varies by 30% (as 
noted in Sec. VIA) even if no further measurements on the 
other loss angle is made. 



B. Torsional Modes of a Coated Hollow Cylinder 

Here we propose an approach with which we can measure 
another combination of loss angles. We consider a cylindrical 
shell with a thin, uniform coating layer outside, as shown in 
Fig. 15(c <^R,d ^ c). In this configuration, the surface defor- 
mations produce strains in the plane of shell according to the 
Donnell shell theory [31]. Here we assumed that there is only 
angular displacement in the shell, which means the longitudi- 
nal position of the cross section won't change. For a torsion 
mode, we only have shear strain energy, the expressions are 
given by 

- (111) 





u 

Us ^dY, (1+g,) 

u cy, (i + (7,) 

As a consequence, the total loss angle can be expressed as 

dY, (1 + C7, 



(112) 



sub " 



(113) 



cF, (1 + (J,) 

For a cylinder shell, according to the Donnell shell the- 
ory, the natural frequency of the n-th torsional mode is given 
by [32] 

1/2 



fn 3 

22L 



P(l 



(114) 



With the measurement of both the thin plate and cylinder 
shell, we can obtain 0b and 0^ of the coating. 



VIII. CONCLUSIONS 

In this paper, by applying the Fluctuation-Dissipation The- 
orem, we obtained a full set of correlation functions (87a)- 
(87c) of Brownian thermal fluctuations of a multi-layer dielec- 
tric coating. In particular, we have related fluctuations of coat- 
ing thickness and coating-substrate interface to independent 
bulk and shear thermal stresses associated with each coating 
layer. These stresses not only induce thickness fluctuations 
of the layers themselves, but also bends the coating-substrate 
interface — this bending noise had not been previously ap- 
preciated intuitively, although its effect has been incorporated 
into formulas, e.g., in Ref. [11]. As a result, we found that 
although thickness fluctuations of different coating layers are 
independent of each other, they each have partial correlations 
with the height fluctuations of the coating-substrate interface. 
Moreover, bulk loss creates a positive correlation between 
them, while shear loss creates a negative correlation. The en- 
tire picture is succinctly written mathematically in Eqs. (89) 
and (90). This coherence structure then gives coating Brow- 
nian noise in Eq. (93). Apart from having provided a peda- 
gogical and systematic derivation of these noise components, 
the most important conceptual consequence of our work is to 
point out an uncertainty in coating loss angles, which has not 
been anticipated previously. We have also incorporated the 
photo elastic effect, the reflectivity fluctuations of the inter- 
faces within the multilayer coating, and considered the effect 
of amplitude modulations caused by Brownian thermal noise. 
All of these turned out to be rather unimportant. 

We have applied our formalism to miiTors that are to be 
used in Advanced LIGO detectors. As estimated in Sec. VI 
and summarized in Table V (calculated for a typical can- 
didate for the Advanced LIGO end test-mass mirror coat- 
ing configuration), parameter uncertainties could lead to non- 
negligible corrections to coating Brownian noise calculations. 
The biggest uncertainties actually arise from the elastic mod- 
uli of coating materials — for example, current uncertainties 
in Young's modulus of the tantala coating material might lead 
up to 60% increase in thermal noise. Although photo elastic 
coefficients for our coating materials are very uncertain, they 
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material 
parameter 


range 


uncertainty 
in 


for details, 
see 


<t>B/<t>S 


0.2 - 5" 


±37% 


Sec. VIA, Figs. 8,9. 


ira 


factor of '-^ 2 


~60% 


Sec. VI B, Fig. 10. 


Cja 


±0.2 


up to 10% if 

te/fe / 1 


Sec. VI B, Fig. 11. 


P 


-KjS <+l 


±1%* 


Sec. VIC, Fig. 12. 



"Fixing the combination (j)o 
'Calculated from Ta205 layers 



Table V: Levels of thermal noise uncertainty caused by parameter 
uncertainties. 



do not significantly affect thermal noise since light does not 
penetrate through many layers. 

It is rather remarkable that our lack of experimental knowl- 
edge about the loss angles, beyond what we had already ob- 
tained from the ring down of drum modes, would not give 
rise to a higher uncertainty in thermal noise. This is rather 
serendipitous, considering our path of understanding of the 
problem: for the baseline parameters of Advanced LIGO, the 
highest contribution to coating Brownian noise arises from 
the coating-substrate bending noise caused by losses in tan- 
tala layers, because these layers are much more lossy than the 
silica layers, and have been assumed to have a much higher 
Young's modulus than the substrate material. This bending 
noise, first elaborated by this work, turns out to be associated 
with the loss angle of the 2-D flexural rigidity, which in turn is 
directly connected to the decay of the drum modes of a thinly 
coated sample. This means the cuiTently existing program has 
been measuring the predominant loss angle all along. Never- 
theless, the level of uncertainty noted in our study still war- 
rants further experiments seeking the other loss angle, e.g., as 
outlined in Sec. VII. In addition, since future gravitational- 
wave detectors may use different substrate and coating mate- 
rials, situations may arise when the loss angle measured now 
does not coiTelate with the total coating brownian noise. 

At this moment, it is worth looking once more at the previ- 
ously used loss angles, 0|| and 0x — although they are mathe- 
matically ill defined, they do correctly reflect the existence of 
two channels of loss. The 0|| was meant to characterize losses 
incurred by the x-y deformations of the coating measurable 
when we do not compress the coating but instead drive its de- 
formations using drum modes of the substrate. This loss angle 
is now replaced by the (mathematically well-defined) imagi- 
nary part of the flexural rigidity, for which extensive measure- 
ments have already been carried out. The was meant to 
characterize the losses incurred by compressing the coating 
layers. This has not been measured because it had not been 
obvious how to easily excite this mode of coating deformation 
(the most obvious way would be to compress the coating layer, 
but that is difficult); however, because the Young's modulus 
of the coating is much larger than that of the substrate, this 
difficult-to-measure loss angle should not contribute as much 
to the total coating noise. This said, in this work, we do come 
up with ways to measure both loss angles, (j)s and (j)B, without 



having to compress the coating layers — but instead by excit- 
ing different modes of substrate deformation. Of course, this 
is only possible because we have assumed that the material is 
isotropic — otherwise we may have to compress the coating 
to directly access the loss induced by such a deformation. 

On the other hand, one may think of the possibilities of 
using substrate materials with higher Young's modulus to re- 
duce the bending noise. Sapphire and Silicon are two viable 
choices because they both have higher Young's modulus than 
tantala. Using Eq. (87a)-(87c), it is straight forward to esti- 
mate the new coating brownian noise while replacing the sub- 
strate material by sapphire or silicon but keeping the same 
aLIGO coating design. It turns out that the coating brown- 
ian noise will be reduced to 35% of its original power spectra 
value if we use silicon substrate or 32% if we use sapphire. 
However, there are other disadvantages for sapphire or silicon 
substrate that prevents us from using them for aLIGO mir- 
rors. The main problem is that they both have very high ther- 
mal conductivity - much higher than fused silica. As a result, 
the substrate thermoelastic noise is one of the important noise 
source for both materials. For instance, if the aLIGO mirror 
was made of sapphire, the bulk thermal elastic noise would 
have about the same magnitude as the coating brownian noise 
at 100 Hz. As for silicon substrate, the bulk thermal elastic 
noise is more than 4 times larger than its corresponding coat- 
ing brownian noise in power because siUcon has even higher 
thermal conductivity than sapphire. One may refer to [37] for 
detailed methods to calculate bulk thermal elastic noise. Set- 
ting up the experiment in a cryogenic enviroment is a possible 
way to reduce the thermooptic noise. 

Furthermore, our formula Eq. (93) can serve as a starting 
point for optimizing the material choice and structure design 
of the multi-layer coating taking light penetration effects into 
account. Our numerical results in Sec. VI D (see Table III) 
have shown that optimization of the coating structure consis- 
tently offers a ^ 6% decrease in thermal noise, regardless of 
In fact, the optimal structure for these ratios are quite 
similar, and configurations obtained for each presumed ratio 
of / 05 are shown to work for other ratios interchangeably. 

Upon completion of this manuscript, we noted that the op- 
timization of coating structure for the case assuming (j)B — (j>s 
(and j3 = 0) has been carried out by Kondratiev, Gurkovsky 
and Gorodetsky [34]. [We note that their formalism is capable 
to treating P ^0 and fc, as well as back-scattering in- 
duced by photo elasticity, but they did not explore the impact 
of these effects in their optimization.] Our results are compat- 
ible with theirs, if we also use these restrictions in parameter 
space and ignore back-scattering. 

A comparison between our result, Kondratiev et al., and 
Harry et al. [11] (which ignores light penetration into the lay- 
ers, and also effectively assumes (j)s = ^b) would therefore il- 
lustrate the effects caused by ignoring photoelasticity and fur- 
ther ignoring light penetration into the coating. This is shown 
in Table VI. This again confirms that for total coating thermal 
noise, light penetration causes noticeable difference in coat- 
ing thermal noise, while photoelasticity causes a negligible 
difference. 
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Coating 


Ref. [11] 
(no light penetration) 


Ref. [34] 
(j6 = and no back scattering) 


This Work 


A/4 


7.18 


7.08 


7.08 


Advanced LIGO 


6.93 


6.82 


6.83 


optimal 


6.73 


6.62 


6.64 



Table VI: Comparison of thermal noise spectral density (assuming (j>B = (j>s and evaluated at 200 Hz, in units of 10 ^'m/\/Hz) between 
different works. 
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Appendix A: Fluctuations of the Complex Reflectivity due to 
Refractive index fluctuations 



Brownian noise is not only caused by the random strains, 
but also by the refractive-index fluctuations caused by such 
strains, through the photo elastic effect [Cf. Eqs. (13) and 
(14)]. We shall quantify this contribution in this section. 



1. The photoelastic effect 

If we denote the displacement of coating mass elements as 
{ux,Uy, M,), then the relative coating-thickness change from its 
equilibrium value can be written as 

51/1 ^u,,, (Al) 

and the relative transverse area expansion can be written as 

dA/A = Mx,x + My.y (A2) 

If we denote 2-dimensional displacement vectors along the x- 
y plane as m {ux,Uy), and two-dimensional gradient as V, 
then we have 



5 A I A = Ux,X + U.yy —V-u 

We can then write the change in refractive index as 



5n = 



dn 


51 


dn 


d\ogl 




(5 log A 



V-M 



(A3) 



(A4) 



where dn/dlogl and dn/dlogA only depends on material 
properties. The two terms on the right-hand side of Eq. (A4) 
represent refractive index change driven by relative length and 
area changes, respectively. The first term is given by [24] 



i3^ = 



dn 
d\ogl 



= -\n^CY 
2 



(A5) 



where C is the photoelastic stress constant, Y is the Young's 
modulus. For silica, CY k, 0.27, therefore )3|[ = —0.41. The 
photoelastic coefficient can also be written as 



n^Pij 



(A6) 



where pij is the photo elastic tensor [36]. Some experiments 
have been done to measure this coefficient for tantala [23]. 
Empirically, the value of pjj varies from —0.15 to 0.45 for 
Ta205 thin film fabricated in different ways. Here for the lon- 
gitudinal photoelasticity, Pj^^ , we use —0.5 in our numerical 
calculation. 

We shall next obtain formulas that will allow us to convert 
fluctuations in n into fluctuations in the complex reflectivity 
of the multi-layer coating. 



2. Fluctuations in an Inflnitesimally thin layer 

Because the coating is much thinner than the beam spot 
size, we only consider phase shifts along the z direction — 
for each value of x. If the refractive index 5n at a particu- 
lar location 5n{z) is driven by longitudinal strain m,, at that 
location, the fact that {uz.z{^)^zz{^')) °^ — z") causes con- 
cern, because this indicates a high variance of 5n at any given 
single point z, with a magnitude which is formally infinity, 
and in reality must be described by additional physics (for ex- 
ample, there would be a scale at which the above-mentioned 
delta function starts to become resolved). Therefore, if we 
naively considers the reflection of light across an interface, at 
z — ZQ, then the independent and high-magnitude fluctuations 
of n{zQ—) and n(zo+) would lead to a dramatic fluctuation in 
the reflectivity 



«(zo-) -»(zo+) 
n(zo+)+n(zo-) 



(A7) 



of the interface, whose magnitude of fluctuation seems to be 
indefinitely large. Fortunately, for any thin layer, if we si- 
multaneously consider propagation through this layer and the 
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Al 









Figure 16: Light propagation across a thin layer (thickness of Al) 
with fluctuating refractive index (from a uniform n2 to an average of 
112 + 5n2 within this thin layer). The propagation matrix correspond- 
ing to this structure is given by Eq. (A8). 



3. The entire coating stack 

Now we are ready to consider the entire multi-layer coat- 
ing. Here we bear in mind that eventually, the fluctuation in 
n has a non-zero coherence length — and we can then divide 
our existing layers further into sub layers with length 51 much 
less than the physical coherence length. Since each of these 
sub layers only makes a negligible contribution to the entire 
complex reflectivity, we only need to consider layers that con- 
tain only one coating material. Let us first focus on Layer 
j, bounded by two interfaces with reflectivities r,_i and rj, 
respectively. The total transfer matrix of the entire stack is 
written as 



reflection and transmission across both of its boundaries, then 
the effect caused by the refractive index fluctuation of this par- 
ticular layer can be dramatically suppressed. Nevertheless, we 
do find an additional fluctuating contribution to the total com- 
plex reflectivity of the multi-layer coating. 

In order to carry out a correct calculation that does not di- 
verge, we first consider a three-layer and two-interface situa- 
tion, as shown in Fig. 16, with ni, «2 and separated by two 
interfaces, with the length of the n2 layer given by Al — and 
here we only consider fluctuations in n2- The entire transfer 
matrix (from below to above, in Fig. 1 6) is given by 



M — Rr|2T02Rr23 



(A8) 



following the same convention as in Sec. II C. Suppose the 
originally uniform «2 now fluctuates, and after averaging over 
this think layer, gives a mean refractive index of n2 + 5n2, we 
use this as the refractive index of the entire layer, and then 
have 



5M-- 



"2 



I —I 

i — / 



5«2 • koM 



(A9) 



This can be considered as a regularization, because each in- 
dividual R,p or R,|, (since their expressions only contain ny, 
nj and n^, but no I) has a standard deviation proportional to 
0(1/a/A/) (when A/ is greater than the coherence length of 
refractive-index fluctuation) or 0(1 /A/) (when AZ is less than 
the coherence length of refractive-index fluctuation) — both 
diverge as AZ — which means the reflectivity fluctuation 
of each of these layers diverge. However, in order for our 
use of average refractive index to make sense in calculating 
the reflectivities r\2 and r23. A/ should be less than the coher- 
ence length of refractive index fluctuations. In any case, the 
total transfer matrix 5M does not diverge; it instead has an 
infinitesimal fluctuation. Moreover, since 5M only depends 
on 5n2 ■ AZ, we shall see that the particular choice of A/ will 
not affect the final results when layers like these are stacked 
together 

The physical meaning of Eq. (A9) is clear: a random field 
of refractive index not only gives a random phase shift (diag- 
onal term), but also gives rise to a random reflectivity (non- 
diagonal term). The latter term is an additional contribution 
that has been ignored by previous calculations. 



M = ---T^^.^,R,T^^R,^._,... 



(AlO) 



Here reflectivity fluctuations within Layer j is going to add to 
the matrix T^^ above. Consider t/z-thick sub-layer located at 
distance z' from the rj boundary (lower boundary in Fig. 1), 
therefore at coordinate location z — Zj+x +z! and integrate, we 
have 



T,|,j T^^. +kQ j 5n {z j+i + z)Tkonjz 



I —I 
i —i 



1 5r]j 
5n* 1 



(All) 



where 



and 



1 f'j 

5nj = - / 5nj{zj+i+z)dz (A12) 

Jo 



577^' = -ikoj 5nj{zj+i+z)e^'''°"j'dz (A13) 
Here we have defined 



TV 

Zj = £ . 

'i=j 



(A14) 



to be the z coordinate of the top surface of Layer j. 

We need to adapt the new transfer matrix into the old form, 
but with modified {rj} and From Eq. (Al 1), since drjj 
is complex, we need to adjust rj, as well as (jij+i'- 



T0,^,Rr,T^^. 



Here we have defined in addition 



(A15) 



drj -tjko J ' 5nj{zj+\ +z) &\n{2kQnjz)dz (A16) 



and 



, r]±l fij 

^Vi = ^ ^0 / dnj{zj+i+z)cos{2konjz)dz (A17) 

•' Irj Jo 
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As we consider photoelastic noise of all the layers together, 
5rj in Eq. (A 16) needs to be used for the effective fluctuation 
in reflectivity of each layer, while 



d = koljdrij + 5x1/ j +S 1/^+ j 



(A18) 



should be used as the total fluctuation in the phase shift of 
each layer. 



4. Unimportance of transverse fluctuations 

Connecting with photoelastic effect, we have explicitly 

5nj{z,x) = Pfujz,x) + j3/V • u (A19) 

Here the vector u is the two-dimensional displacement vec- 
tor {ux-Uy) and V- is the 2-D divergence along the x-y plane. 
For terms that contain u, we note that when taking the optical 
mode into account [see Sec. II D], i.e., when a weighted aver- 
age of ^ is taken, they yield the following type of contribution 



j Iix)(v-u 
dl (n • ul) 



cf'x 



dM 



u- VI d X 



M 



M 



(A20) 



Here M stands for the 2-d region occupied by the beam, and 
dM is the boundary on which power already vanishes. As a 
consequence, the first term is zero according to the boundary 
condition, while the second term gains a factor of (/,/rbeam) 
with respect to other types of coating Brownian noise, here 
Ij is the thickness of the j-th layer, and rbeam is an effective 
beam radius. Since we always assume coating thickness Z, to 
be much smaller than the beam radius rbeam; we can neglect 
refractive index fluctuation due to area fluctuation. 



Appendix B: Elastic Deformations In The Coating 

Throughout this paper, we assume the mirror substrate to 
be a half infinite space. We establish a Cartesian coordinate 
system, with {x,y) directions along the coating-substrate in- 
terface, and z direction orthogonal to the mirror surface (in 
the elasticity problem, we also ignore mirror curvature). This 
allows us to calculate elastic deformations in the spatial fre- 
quency domain. We will also assume the coating thickness to 
be much less than the beam spot size. 

We denote the displacement along x, y and z directions as 
Ux, Uy and u^. It is then straightforward to express the 3x3 
strain tensor S in terms of their derivatives, and stress tensor 
T in terms of Hooke's Law: 







1 / dui 

2 \dxj 

S.i 
1 



dxj 



[S.j+Sji] 



-5,e 



(Bl) 
(B2) 
(B3) 
(B4) 




Figure 17: Sample with single layer coating, force is applied perpen- 
dicular to the air/coatine interface. 



Here we have x^ — {x^y^z), with Latin indices (like / and /') 
running from 1 to 3. Within any layer, it is straightforward to 
write down the most general solution of the elasticity equilib- 
rium equation 







= ikx[(a+ + Kz^+)e^~ + (a_ - Kz^- )e 

- iky[y+e''' + y_e-'''] 

= iky[{a+ + Kzp+)e'^' + (a_ - Kzp-)e' 

+ ikxif+e^' + f-e-"-] 

= -K-[a++j3+(-3 + 4(7 + fcz)]e''^ 

+ K-[o:_+j3_(-3+4(7-K-z)]e"''^ 



(B5) 

(B6) 
(B7) 
(B8) 



where tilde denotes quantities in the x-y spatial-frequency do- 
main, and K — ^/k^+k^,. Namely 



(B9) 



We now consider a single-layer coating on a substrate, 
with the coating-substrate interface located at z 0, and the 
coating-air interface at z = /. Suppose there is a force profile 
F{x,y) exerted perpendicular to the coating surface, at z = li, 
< d < I, and let us calculate the elastic deformation field 
caused by F. The entire system is now divided into three re- 
gions, (a): d <z< I, (b): < z < d, and (s): z < 0. At the 
interfaces, we obtain the following 15 boundary conditions. 



'T'CI 

iz 


= 0, 


z 


= 1 


(BIO) 






z 


= d 


(BID 






z 


= d 


(BI2) 






z 


= 


(BI3) 
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as well as the condition that when z ^ (which 

leads to = = fL ~ 0). We are left with 15 fields of 

(a^J^,g,aiJi,}^,a|,j3|,n) (B14) 



Using linear superposition, as well as taking the appropriate 
limits of the above solution, it is straightforward to obtain 
elastic deformations in all the setups required in order to ob- 
tain cross spectra between different noises. 



which can be solved from the 15 boundary conditions. As- 
suming Kd I and fc/ ^ 1, we obtain that all 7 vanish, and 



PI 



at 



Pi 



F{1 + g,)[2 - 3g, + g,(-3 + 4g,)] 
2F,)c2(-l + (7,) 

F(g,-g,)(l + g,) 

2F,K-2(-l + (7,) 

F(l + g,)(-3+4g,) 

4y,fc2(-i + (j,) 



4y,K-2(l-(J,.) 

F{1 + g,) [2 - 3g, + gc(-3 + 4g,); 
2y,K-2(-l + (7,) 

F(g,-g,)(l+g.v) 

+ g) - y,(-3 + gy +4g2)] 

4FF,K-2(-l + (7,.) 

F[y,(i + g,)-y,(i + g,)] 
4yy,K-2(-i + cj^) 

F(l + cj,)(-l+2a,) 



(B15) 
(B16) 
(B17) 
(B18) 
(B19) 
(B20) 
(B21) 
(B22) 
(B23) 
(B24) 



We can therefore obtain the stain tensor in the frequency do- 
main for the coating, the non-zero elements for region (a) are 
given by 



F^,2(_i+2g,)(l + g2) 

F,K-2 



Fk,ky{-l+2a,){l + a^) 



YsK^ 

^ g,(-l + g, + 2g2) 

y.(-l + CT,) 

while those in region (b) are given by 

Ffe2(-1+2(J,)(1 + CT2) 



SI 
St 



Sxy — Syx 



F^2(_i+2(J,)(l+g2) 

y,K-2 

^ Fk,ky{-l+2a,){\ + a}) 

Y,K^ 

-(l+2g,) g,(-l+(y,+2g2) 

Yc y.(i-c7,) 



(B25) 
(B26) 
(B27) 
(B28) 

(B29) 
(B30) 
(B31) 
(B32) 



Appendix C: Definition of loss angle 

In the past [11], the coating loss angle was defined in as- 
sociation with the parallel and perpendicular coating strains. 
The equation is written as 



coated 



'sub 



-01 



(CD 



U u 

where 5U\\ and 8U± are the dissipation energy density in par- 
allel and perpendicular coating strains 

5f/|| = ^ {S„T„ + SyyTyy) dxdj (C2) 

5t/± = ^ ]^S-,J,-, dxdy (C3) 

and where Sij are the strains and 7]y are the stresses. While 
such a definition seems to be compatible with the symmetry 
of the system, the quantities 8U\\ and 5f/x cannot be used 
as energy, since in certain scenarios they each can become 
negative. 

For example, if we have a cube with surface area of each 
side A (poisson ratio a. Young's modulus Y), and we uni- 
formly apply two pairs of forces, one pair with magnitude / 
on opposite yz planes, the other with magnitude F on oppo- 
site xy planes, with f <^F, as shown in Figure 18. According 
to definition of Young's modulus and Poisson's ratio, up to 
leading order in f/F the non-vanishing strains are. 



5zz = 



F/A 



SxX Syy CT 



F/A 



(C4) 



On the other hand, for stress, we have, up to leading order in 

f/F, 



-f/A, Tyy^O, 



-F/A. 



As a consequence, we have 



= S,,T,, + SyyTyy = ~afF/{A^Y) < 



(C5) 



(C6) 



which means 8U\\ is not a reasonable candidate for energy, at 
least with a y^O. Since it is also true that S^xTxx < we will 
arrive at 



(C7) 



if we take this configuration and rotate for 90 degrees around 
the y axis, such that x rotates into z- 

One reasonable way of defining the loss angle is to derive 
from the fundamental elastic energy equation. The general 
form of the stored elastic energy density U can be written as 



Us = 



(C8) 

(C9) 
(CIO) 
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iF Appendix D: Advanced LIGO style coating 




Figure 18: Solid cube witli two pairs of forces applied on the side: 

Where K is called the bulk modulus and /i is the shear modu- 
lus. In the calculation, we use Young's modulus Y and Pois- 
son's ratio a instead of K and /i. Their relation are given in 
Eq.(15). The expansion © and shear E are both irreducible 
tensorial parts of the strain tensor S. 

® = Su (Cll) 
I.= ^-[S,j+Sj,)-^-gijSkk (C12) 

Note that the expansion and shear energy Ug and Us is always 
positive, so it is consistent to define loss angle by 0b and ^s- 



j 



1-5 


0.497325 0.208175 0.289623 0.237274 0.250176 


6-10 


0.245330 0.249806 0.240129 0.270968 0.224129 


11-15 


0.251081 0.259888 0.260826 0.213460 0.290468 


16-20 


0.214524 0.273240 0.230905 0.259924 0.230020 


21-25 


0.275429 0.233086 0.270385 0.208581 0.273798 


26-30 


0.249741 0.267864 0.204698 0.292317 0.209712 


31-35 


0.278560 0.220264 0.282694 0.221687 0.268559 


36-38 


0.233460 0.270419 0.223050 



Table VII: Structure of an Advanced LIGO-like coating optimized 
jointly for dichroic operation and thermal noise. Thickness of 
each layer is given in units of wavelength (for light with vacuum- 
wavelength of 1064 nm) are listed here for the 38 layers. Note that 
/l 3 5 are Si02 layers, while hAfi.... are Ta205 layers. 



In Table VII, we provide the structure of the coating opti- 
mized jointly for dichroic operation and thermal noise (base- 
line parameter, neglecting penetration). 
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